﻿ Lucrarea nr 6 Transmisii de date în banda de bază cu interferenţă controlată intersimbol Scopul lucrării: se realizează simularea trecerii unei secvenţe de date bipolare printr- un canal de comunicaţii utilizând ca filtru global un filtru de tip “cosinus” Se simulează astfel o transmisiune cu interferenţă controlată inter-simbol 6 1 Transmisii de date în banda de bază cu interferenţă controlată intersimbol Avem în vedere sistemul de transmisiuni în banda de bază, din figura 1: EMIŢĂTOR CANAL RECEPTOR z(t) a(t) s(t) s’(t) r(t) {âk} {ak} Eşantionare/ G(ω ) MIA E) GR(ω C(ω + ) Decizie tk Figura 1 Sistem de transmisiuni de date în banda de bază unde : GG CGωωωω=⋅⋅ (1) () ()()() ER reprezintă răspunsul în frecvenţă al lanţului de transmisiuni emiţător-canal-receptor Cω=1 (cu fază (caracteristica globală a sistemului) Canalul se consideră ideal, () nulă), ce nu perturbă semnalul transmis prin el În cazuri reale însă, canalul introduce o serie de distorsiuni asupra semnalului: întârzieri, atenuări, zgomote În cele ce urmează vom negija efectele de întârziere şi de atenuare ale canalului, ţinând cont doar de zgomotul introdus de el, pe care-l vom considera ca fiind unul aditiv, alb şi cu ) 1=ωC distribuţie gaussiană, cu dispersia constantă σ i Astfel, vom considera ( Influenţa zgomotelor introduse de canal se pune în evidenţă prin însumarea ) ts, cu un “zgomot”, z(t) semnalului de la ieşirea canal, ( ωG asigură rata de Modelul “filtru trece-jos ideal” al caracteristicii globale () transfer maximă a transmisiei datelor fără interferenţă intersimbol, în condiţia unei limitări a benzii de frecvenţă (teorema lui Nyquist) Însă modelul este unul ideal ce nu corespunde situaţiilor reale Celelate modele stabilite, familia de filtre tip “cosinus- 1ridicat”, permit, la aceeaşi bandă de frecvenţe disponibilă, o rată de transfer maximă mai redusă Se poate arată că, permiţând un anumit grad de interferenţă între simbolurile transmise, se găsesc şi alte caracteristici de transfer ce permit rata maximă permisă de teorema lui Nyquist în condiţii similare de bandă Un astfel de model este cunoscut sub denumirea de “filtru global tip cosinus” Ră spunsul în frecvenţă al sistemului este dat de relaţia: T2 πω ⎧ cosT ≤ω⋅ ⎪ T4 () G (2) =ω ⎨ π 2 ⎪ >ω 0 ⎩ T 2 limitat în bandă la pulsaţia π T Răspunsul la impuls al filtrului este: t 2π cos ∞ 41 ω tj T () () deGtg (3) ⋅=ωω= 2∫ ππ t 42⎞ ∞− ⎛ 1⎟ − ⎜ T ⎠⎝ ) ωG, din ec (1), cu În figura 2 a se prezintă funcţia răspuns în frecvenţă, ( valoarea T=2, în timp ce în figura 2 b se prezintă funcţia răspuns la impuls, g(t), corespunzător (generat în N=128 puncte) a) b) Figura 2 a) Răspunsul în frecvenţă al filtrului “cosinus”; b) Răspunsul la impuls corespunzător Funcţia răspuns la impuls g(t) verifică: ()() 144TgTg= −= (4) () 10024,kTkTg≠∀=⋅+− Această proprietate permite o transmisie de date cu durata 2TT= 0 a simbolurilor Astfel, dacă se consideră un semnal a(t), de intrare,ca fiind un mesaj ce constă din N a, ce apar la intervale de timp T, de forma: simboluri consecutive, {} k0 2 −1N () −δ=nTta)t(a (5) ∑ 0n =0n ) ωC ce nu introduce zgomot atunci, semnalul recepţionat, r(t), în cazul canalului ( peste semnalul util, are expresia: −1N () −⋅=nTtga)t(r (6) ∑ 0n =0n 0T , avem: Eşantionând semnalul r(t) în momentele de timp kTt+ = 0k− 2 1−N TTT ⎛⎞⎛⎞⎛ 000⎞ () ++−=⋅−−−⋅== ⎟⎜⎟⎜⎟⎜ ∑ 10− kknkkgagaTkngar)t(r ⎠⎝⎠⎝⎠⎝ 0222=n (7) 1−N T ⎛ 0⎞ () +=⋅−−−⋅+aaTknga ⎜ ⎟ ∑ 10−kkn ⎝ ⎠ 102−≠=k,kn,n deci eşantionul prelevat la un moment t, r, depinde atât de simbolul curent cât ş kki de cel anterior transmis 6 1 1 Influenţa zgomotului asupra proiectării sistemului În toată descrierea anterioară, ignorând zgomotul, ne referim doar la răspunsul la impuls al ansamblului filtru de emisie-canal-filtru de recepţie, fără a insista asupra formei particulare a răspunsurilor la impuls ale celor două filtre (emisie şi recepţie) În cele ce urmează se va arăta că prezenţa zgomotului în canal poate fi combătută printr-o anumită modalitate de alegere a acestor filtre tg răspunsul la impuls al filtrului de emisie, la ieşirea acestuia se Notând cu () E obţine semnalul: −1N −⋅=nTtga)t(s (8) () ∑ 0En =0n semnal ce este transmis prin canal Canalul fiind idealizat, deci fiind unul ce nu atenuează şi nici nu întărzie sau dispersează, la intrarea filtrului de recepţie vom avea semnalul: ()()() tztst's+= (9) relaţie ce pune în evidenţă cele două componente ale unui semnal ce trece i componenta de zgomot, printr-un canal zgomotos: componenta de semnal util, s(t), ş z(t) Filtrul de recepţie fiind un sistem liniar rezultă că şi semnalul de la ieşirea acestuia se poate descompune în două componente, una de semnal util şi cea de zgomot Astfel, semnalul recepţionat, r(t), se poate scrie: 3 )()() ( tztrtr+= (10) u0 spunsul filtrului de recepţie unde () tr reprezintă u componenta de semnal util, deci ră tz este componenta de zgomot, deci răspunsul filtrului de la semnalul s(t), iar () 0 recepţie la zgomotul z(t) din canal Componenta de semnal util are expresia: −1N () −⋅=nTtga)t(r (11) ∑ 0nu =0n 0T Eşantionând semnalul r(t) în momentele de timp kTt+ , avem: = 0k− 2 () () +==+= kukkkkukzrrtztr)t(r=0 N−1 TTT⎞ ⎛⎞⎛⎞⎛ 000 () +−=+⋅−−−⋅= ⎟⎜⎟⎜⎟⎜ ∑ kkkngagazTknga+− 10 (12) ⎝⎠⎝⎠⎝ n⎠=0222 N−1 T⎞ ⎛ 0 () +=+ ⋅−−−⋅+ ⎟ ⎜ ∑ kkkknzaazTknga+−10 ⎝ k,kn,n⎠−≠=102 relaţie identică cu (4), cu excepţia zgomotului La recepţie, decizia asupra simbolului binar transmis, a, notat aˆ, se poate kk face pe baza eşantionului curent r ş condiţionare a ki deciziei anterioare 1−kaˆ Această deciziei curente de decizia anterioară face posibilă “propagarea” unei decizii incorecte (ce se manifestă prin luări de decizii incorecte până la apariţia unei erori noi) În practică această posibilitate se elimină prin precodarea secvenţei de date ce ase transmite secvenţ urmează a fi transmis Astfel, la emisie, în loc de secvenţa { } ka {} de: kbdată ⊕ =bab (13) 1−kkk Se defineşte raportul semnal/zgomot, RSZ, la ieşirea filtrului de recepţie la momentul t ca fiind: 0 2 () 0utr RSZ = 0 (14) P 0z irea filtrului de recepţie unde P este puterea zgomotului de la ieş 0z ) Se caută filtrul de recepţie optim, ( ωG, care să maximizeze raportul R semnal/zgomot la ieşirea filtrului de recepţie la momentul 0= 0t Acest filtru se găseşte ca fiind unul cu expresia răspunsului în frecvenţă: 4 * )() (ω⋅=ωSkG (15) R deci, filtrul adaptat la forma semnalului din canal Necunoscând însă, apriori, informaţia ce se va transmite prin canal, cerinţa (14) este dificil de realizat O soluţie mai puţin pretenţiosă, insă doar suboptimală, acceptată în practică ca şi una de compromis, este oferită de condiţia de maximizare a raportului semnal/zgomot la ieşirea filtrului de recepţie în cazul transmiterii doar a unui singur simbol de informaţie Acest lucru duce la o expresie forma: 2 ()()() ω⋅=ω⇒ω⋅=ωGkGGkG (16) () EER Condiţie satisfăcută, cu k=1, de alegerea: ()() ω=ω=ωGGG (17) () ER Relaţia (14) nu impune nici o restricţie deosebită caracteristicii de fază a celor două )() funcţii răspuns în frecvenţă, ( ωG respectiv ωG, singura condiţie ce trebuie să ER fie verificată este acea de sumă nulă a celor două faze la orice frecvenţă Pentru o funcţie () ωG reală, cum este şi în cazul nostru, putem alege: ()() ω=ω=ωGGG (18) () RE Într-un sistem binar monopolar de transmisiuni, simbolul a poate fi 0 sau 1, în timp k ce un sistem binar polar va emite simboluri a cu nivelele posibile -1 ş ki 1 Valorile in tabelul 1 posibile ale eşantioanelor r, în cele două k cazuri, se prezintă Tabelul 1 Valorile posibile ale eşantioanelor r, în funcţie de simbolurile transmise k Sistem unipolar Sistem bipolar ka 1− ka kr ka 1−ka kr 0 0 0 -1 -1 -2 +1 0 +1 +1 -1 0 0 +1 +1 -1 +1 0 +1 +1 +2 +1 +1 +2 La recepţie, decizia asupra simbolului binar transmis, a, notată k kaˆ, se poate ie face pe baza eşantionului curent r ş ki a deciziei anterioare 1−kaˆ Metoda de detecţ utilizată permite deasemenea detectarea unor decizii incorecte Astfel, de exemplu, în cazul unei transmisii bipolare, o valoare r de -2 nu poate fi urmată k de o valoare facilitate însă nu elimină posibilitatea luării unor decizii incorecte, 2+ 1=+kr Această mai mult, din cauza condiţionării deciziei curente de decizia anterioară, apare posibilitatea “propagării” unei decizii incorecte (luări de decizii incorecte pănă la apariţia unei noi erori de decizie sau apariţia unei violări de semn) ului de transmisiuni în banda de bază 6 2 Procedura de simulare a lanţ 5 Echivalând sistemul de transmisiuni în timp continuu din figura 1 cu un sistem în timp discret şi discretizând, deasemenea, semnalele de intrare/ieşire, se obţine structura sistemului de analizat din figură: EMIŢĂTOR CANAL RECEPTOR ad[n] Sd[n] sd’[n] rd[n] {âk} Eş {ak} antionare/ C(Ω) ) GR(Ω MIA GE(Ω ) Decizie nk Figura 3 Sistemul în timp discret echivalent sistemului din figura 1 Procedura de simulare implică: 2 1 Discretizarea semnalului de intrare (at, prin eşantionare cu pasul T Se ) e obţine astfel o secvenţă numerică [ an anT=, unde se alege () ] de l ,l,LTT02>==, rezultând o lungime totală MNL=⋅ a secvenţei Alegând be0b a lui 2, va rezulta o numărul de simboluri de transmis N ca fiind o putere întreagă a lui 2 lungime totală a secvenţei, M, tot o putere întreagă 2 2 Echivalarea sistemului în timp continuu cu un sistem în timp discret pe baza metodei de echivalare prin invarianţa răspunsului la impuls Se găseşte funcţia GΩ, funcţie ce, în condiţiile răspuns în frecvenţă a sistemului numeric echivalent, ( ) d l 02>==l,LTT (ce respectă teorema eşantionării), va avea expresia: 0be Ω π ⎧ bL Ω=≤ΩcosL ⋅2 ⎪ 0b () ( ) π+Ω=ΩGG,LG (19) =Ω22 () ⎨ ddbd ⎪ π<Ω<Ω 0 0 ⎩ În vederea reprezentării pe calculator a acestei funcţii ea se discretizează, 2 ⋅=Ω, rezultând o secveţă numerică eşantionând uniform în punctele kπ Mk Gk G=Ω , cu perioada M (o putere întreagă a lui 2) Secvenţa se [] () ddk periodică memorează reţinând M valori succesive, (o perioadă), ale ei seşte secvenţa răspuns, [rn, de la ieşirea filtrului de recepţie 2 3 Se gă ] d Ţinând cont de faptul că, prin ipoteză, canalul nu introduce zgomot peste semnalul util, aceasta se calculează în domeniul frecvenţă, prin: - găsirea secvenţei [] A k FFT a n=, calculată în M puncte; []() dd - calculul produsului [] Rk Ak Gk=⋅, tot în M puncte; [][] ddd - calculul secvenţei [] rn IFFTRk= [] () dd 6 bL 2 4 Găsirea eşantioanelor [rrn=, cu Lkn⋅ , ce corespund −= ] kdkbk+ 2 rt, de la ieşirea filtrului de recepţie din eşantionării semnalului în timp continuu,( ) 0T = Eşantionarea se poate face înmulţind figura 1, în momentele de timp kTt+ 0k− 2 rn, cu o secvenţă de “eşantionare”, [pn , de forma: secvenţa [] d] () + LLn21M ⎧ bb [] =np (20) ⎨ restin0 ⎩ 2 5 Luarea deciziilor asupra simbolurilor a transmise pe baza eşantioane-lor k r În acest scop se scade, din eş tor decizei kantionul curent rk, nivelul corespunză anterioare, aˆ, iar rezultatul, 'r , se compară 1−kk cu un prag d (1/2 în cazul unei secvenţe unipolare şi 0 în cazul unei secvenţe bipolare), şi se decide simbolul: ≥ d'rdaca1 ⎧ k = aˆ (21) ⎨ k < d'rdaca0 k ⎩ 3 Partea experimentală 3 1 Utilizând un program de calcul ş tiinţific (Matlab) găsiţi secvenţa răspuns Gk, în M=128 puncte cu în frecvenţă a sistemului în timp discret echivalent, [ ] d 16=L b siţi şi reprezentaţi funcţia 3 2 Aplicând algoritmul FFT invers, (IFFT), gă răspuns la impuls, [] ng, al acestui sistem Găsiţi poziţiile în care secvenţa ia d valoarea 1 Cât rezultă intervalul dintre aceste poziţii ? ie de generare de valori aleatoare, generaţi o secvenţă 3 3 Utilizând o funcţ ] na, de lungime M=128, versiunea discretizată a unui aleatoare de date unipolare, [ d semnal a(t) de 8= =LMN simboluri unipolare b rind procedura descrisă în §2 3, găsiţi secvenţa []nr ce se obţine la 3 4 Urmă d ieşirea sistemului numeric dacă la intrarea sa se aduce secvenţa [] na În ce domeniu d ] de valori se încadrează valorile [ nr ? Justificaţi d i momentele de timp optime de eşantionare, n, şi găsiţi secvenţa de 3 5 Stabiliţ k $ biţi obţinuţi în urma deciziei, {} a, conform procedurii descrise în §2 4 ş ki §2 5 Ce se poate constata ? Ce valori au eşantioanele la momentele Lkn⋅ ? bk= i 3 3-3 5 pentru cazul unei secvenţe de date bipolare 3 6 Repetaţ 7